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ABSTRACT 


A theoretical analysis was conducted of the dynamic behavior of micron 
size particles entrained in gas flow on the two-dimensional blade-to- 
blade surface of a circular stationary cascade of turbine stator blades. 
The particle velocity lag and angular deviation relative to the gas 
was determined as a function of particle diameter and mass density. 
Particle size and density were varied over ranges selected to corres- 
pond to typical laser-Doppler velocimeter (LDV) flow field mapping 
applications. It was found that velocity lag and angular deviation 
increased whenever particle size or mass density increased, and that 
particle tracking was more sensitive to a change in particle diameter 
than to a change in mass density. Results indicated that LDV applica- 
tions employing 1 gm/cc tracer particles with diameters greater than 
approximately 1 micron, or 0.5 micron diameter particles with mass 
densities greater than 4 gm/cc would experience velocity and angular 
deviations generally greater than 2 percent and 1 degree, respectively. 
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s trea.ni— channel thickness normal to meridional streamline 
drag coefficient 

mass flow rate of gas through stream channel 
Mach number 
particle radius 

radius from axis of rotation to meridional stream-channel mean line 

Reynolds number, 2rp ju-V|]i 

6 

temperature 

gas velocity 

particle velocity 

meridional streamline distance 
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gas viscosity 
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absolute angular deviation between the particles and gas, I 3p $g[ 
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INTRODUCTION 

Laser-Doppler velociinetry (LDV) is a noncontact, optical tech- 
nique for determining the velocity of flowing fluids by measuring 
the velocity of tracer particles entrained in the fluid rather than 
by measuring the velocity of the fluid itself. The LDV technique 
is of considerable importance in gas dynamic measurements because it 
is unnecessary to introduce any hardware into the flow system. 

Because of the absence of a physical probe, the LDV technique is ideally 
suited to measuring gas velocities in the flow passages between rotor 
and stator blades of turbomachines. Development of noncontact vel- 
ocimetry techniques in turbomachinery will provide data currently 
required to extend and improve design methods blade performance, 
and efficiency of components used in advanced airbreathing engines. 

It is for this reason that the LDV technique is beginning to find ap- 

1 O 

plication in the aircraft turboTnachinery field. 

The laser-Doppler velocimeter system measures only the motion 
of micron size naturally occurring or artificially generated tracer 
particles entrained in the fluid, and its accuracy is, therefore, 
limited by the accuracy with which the particles follow the fluid flow. 
Since it is important to know how accurately the particles track the 
fluid flow, careful consideration must be given throughout the flow 
region to the magnitude and direction of the particles velocity rel- 
ative to the carrier fluid velocity. This paper documents a theoretical 
investigation of the dynamic behavior of micron size particles entrained 
in gas flow on the two-dimensional blade-to-blade surface of a circular 
stationary cascade of turbine stator blades. The particle velocity 
lag and angular deviation relative to the gas are determined as a 
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function of particle size and mass density. Information of this nature 
is currently required for the application of LDV techniques to the 

intrablade flow field mapping of aircraft turbomachinery. 

3-5 

Several previous studies have investigated particle tracking 
for LDV applications but have been restricted to particle flow in. such 
systems as turbulent flow's, nozzles, oscillating gases, shock waves, 
and over airfoils. Reference 6 used a LDV system to measure the 
velocity lag of aerosol particles in small supersonic nozzles and 
then used this information in conjunction with numerical predictions 

7 

to determine particle size. A theoretical study by Hussein and Tabakoff 
and an experimental study of Tabakoff et al,® provide trajectories and 
velocities of solid particles suspended in a gas passing through an 
axial flow turbine stage. The objective of their investigations was 
to estimate the flow properties of a gas-particle suspension flowing 
over turbine blades. This information would ultimately be used to 
minimize blade erosion, through appropriate redesign. Since their work 
considered particle sizes generally ranging from 100-900 microns, their 
results are inapplicable for LDV applications because of the resulting 

r\ 

high velocity lags. Wisler and Mossey^ obtained detailed flow field 
measurements within the rotating blade row of a low speed axial com- 
pressor using a laser-Doppler velocimeter. Their measurements generally 
agreed to within 2-3 percent with hot film and pitot-static probe 
measurements outside the blade row, and with potential flow solutions 
near the leading edge. Their work was restricted however , to flow 
velocities sufficiently low that particle velocity lag was not regarded 
as a significant problem. However, as gas flow velocities and accel- 
erations increase, particle velocity lag in blade channels becomes 
more significant and application of LDV techniques to flow field analysis 
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should be guided by the results of a dynamic analysis of similar 
gas-particle flow systems. 

GAS-PARTICLE FLOW SYSTEM 

The flow system analyzed in this study is the gas -particle flow 
on the mean blade-to-blade surface of a circular stationary cascade 
of turbine stator blades. The cascade was designed for cooled-tur- 
bine thermodynamic studies and is currently undergoing testing at the 
NASA Lewis Research Center. It has 72 blades and each blade has a 
3.85 in. span and 2.47 in. chord. The tip diameter of the turbine is 
31,9 in.; the hub-to-tip radius ratio is .76 and the mean radius 
solidity is 2.05. Coorinates of the actual blade profile are given 
in ref. 9 for three radial sections. The flow Is considered to be 
two-dimensional and the independent variables are the meridional 
streamline distance X and the tangential angle 0 . For purposes of 
analysis, a finite solution region is chosen on the blade-to-blade 
surface. It is bounded on the top and bottom by adjacent blade sur- 
faces , and by upstream and downstream boundaries, along which the flow 
is assumed to be uniform.-. The blade-to-blade sumface of revolution, 
solution region, and coordinate system employed are shown in fig.l. 

A stream channel is defined by specifying a meridional streamline 
radius R arid a stream- channel thickness b at several meridional locations 
The streamline radius and stream-channel thickness are R=l. 16773 feet 
and b=. 0031939 feet, respectively, and remained constant (to within 
1 percent) throughout the solution region. The Ib'lade -to -blade solution 
region and the grid used for finite -difference amalysis of the governing 
equations is shown in figure 2. 
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GOVERNING EQUATIONS AND ANALYSIS 

Analysis of the gas-particle flow field was conducted in two 
steps. In the first, the magnitude and direction of the gas velocity 
was determined at all Interior grid points of the blade-to-blade 
solution region. In the second step, the magnitude and direction of 
the particle velocity was determined throughout the region by analyzing 
the path followed by particles in a two-dimensional flow field with a 
known gas velocity distribution. Decoupling the gas and particle flows 
in this way is justified because of the low mass fraction of naturally 
occurring or artificially generated aerosols In most LDV gas-dynamic 
applications. Therefore, the particles do not influence the gas flow, 
and the particles and gas may be analyzed separately. 

GAS FLOW 

The gas flow field in the blade-to-blade solution region was 
determined with a solution method developed by Katsanis 10 . The method 
yields a two-dimensional , isentropic, shock free, transonic flow so- 
lution on a blade-to-blade surface on which the flow is essentially 
subsonic. The solution is obtained by a combination of a finite- 
difference stream-function solution and a velocity gradient solution. 
The simplifying assumptions used in formulating the governing equation 
are: (a) the fluid is a perfect gas and is nonviscous , (b) the 

flow is steady and irrotational (c) the velocity component normal 
to the solution surface is zero, (d) the magnitude and direction 
of the velocity are uniform across the upstream and downstream 
boundaries, (e) the only forces acting on the fluid are those due to 
momentum and the pressure gradient, and (f) the total energy of 
the system remains constant. 
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The conservation of momentum equation for a compressible, non- 
viscous fluid in steady flow through a stationary cascade is given^^ 
in terms of the stream function as: 

ijlS? + -L I 'kpii'c/y ■ bM a’mo C i_ ( b phi s 6 (l) 

bX 1 R 4 oS" R*p a ii b® bx L R bp s bX ' 3 i 


where 




K b ^ b© 

U 0 - ^ 

0 bp 3 bX 


The stream function is 0 on the upper surface of the lower blade 
(suction surface) and 1 on the lower surface of the upper blade 
(pressure surface). For flow systems in which the flow is everywhere 
subsonic, the solution to eqn. (1) can be obtained with a finite- 
difference approach. However, if there is locally supersonic’ 
flow, the solution is obtained by a velocity-gradient method using 
information obtained from a finite-difference solution of eqn. (1) at 
a reduced weight flow. The derivation of the velocity-gradient equa- 
tion and details of the solution precedure are given in ref 10. 


PARTICLE FLOW 

The differential equation of motion in vector form of a single, 
spherical particle immersed in a fluid in the absence of potential 
force fields is 


<t) 





The term on the left denotes the mass times the acceleration of 
the particle, and, therefore represents the force required to accel- 
erate the particle. The term on the right represents the total viscous 
drag exerted on the particle by the surrounding fluid and is ex- 
pressed in terms of an empirically determined viscous drag ceofficient 
Cp. Equation (2) assumes that the force on the particle caused by 
the pressure gradient in the surrounding fluid, the force required 
to accelerate the apparent mass of the particle relative to the fluid, 

and the force which accounts for the deviation of the flow from steady 
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state (Basset force) are negligible. These assumptions are . jus tified 
in flow systems whose properties are changing very gradually or when 
the particle density is much larger than the fluid density. Both of 
these conditions are satisfied in this study. It is also assumed that 
the particles are non-interacting and uniformally distributed within 
the carrier gas, and that there are sufficiently few particles that 
their presence does not alter the gas properties from that of the gas 
flow alone. These assumptions are valid for most LDV applications 
employing naturally occurring or artifically generated aerosols. The 
steady-state meridional and tangential components of eqn. (2) 
become 


dV* = 
dx 

( 0 X lH-vl 

BrPpV x 

(3) 

dV e _ 
de 

(Ue-V^IU-Vl 

8rp p V 0 

(A) 


The path followed by a particle and its velocity along that path 
may be determined by numerically integrating eqns 0 (3) and (4) in 
an interative manner throughout the gas flow field. 
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Numerical integration is required because of the nonlinear nature of 
the empirically determined drag coefficient. Since the gas flow is 
assumed not to have a radial velocity component, the particles also do 
not have one, and are, therefore, constrained to remain in the stream 
channel defined by the meridional radius R and the stream-channel thick- 
ness b. 

In the absence of reliable drag data the drag coefficient is fre- 
quently assumed to be given by the classical Stokes drag law “ 24/Re, 
where Re = 2rp |u-V|/p. Although It is most convenient to use Stokes 
law for the drag coefficient, it becomes increasingly unrealistic for 
Reynolds numbers above approximately 1.0 and in flows where rarefaction, 

compressibility and inertial effects are significant. The drag coe- 

11 

fficient used in this study is given by Carlson and Hoglund and corrects 
the Stokes law for these effects. Therefore, it is applicable throughout 
the continuum, slip and transition flow regimes. The expression is 


r , Ml , , , 4.lsS # 1 

2^ 1 (l4»»i%Re ) ( 1 4- -axp (- /frA - 3/Pa. )") 

ReL i + (wv/rO( 3 + l.iS exp C-i.isRe/M)) J 


where the term in the bracket accounts for deviation from Stokes flow. 
The Mach number is based upon the absolute relative velocity of the 
gas with respect to the particles. 

PARAMETRIC STUDY 

The dynamic behavior of several gas -particle mixtures was deter- 
mined by numerically solving eqns . (3) and (4) in an interative manner 

for the particle velocity components. The viscous drag coefficient is 
given by eqn (5), Several cases were examined by varying the size and 
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mass density of the particles over ranges selected to correspond to 
values encountered in most LDV gas dynamic applications. The par- 
ticle radius varied from 0.5 - 2 microns, and the mass density ranged 
from 1-4 gm/cc. The gas flow remained constant throughout the study 
and was characterized by a mass flow rate, rn = 6.68 x 10 3 lb/sec, 
mass density, p = 7.403 x 10 -2 lb/ft 3 , inlet temperature equal to 

g 

523,7 R, and a specific heat ratio y=l o 40. The gas flow entered the 
solution region with an inlet flow angle of zero relative to the meri- 
dional direction and an average velocity of 285.4 ft /sec. Particles 
were introduced into the solution region at several locations along 
the upstream boundary. They were. assumed to possess a velocity equal 
to 99 percent of the corresponding gas velocity at that point, and have 
no angular deviation relative to the gas. Their trajectories and 
velocities were numerically determined as they progressed through the 
region. Analysis of a particular pa rticle was continued until either 
it passed out of the solution region or until it entered the near-blade 
region. This region is defined as the area between the finite-differ- 
ence grid network and a blade surface. Since particles which entered 
this region were generally large and dense, and therefore unsuitable 
for LDV applications, their analysis was terminated ^immediately after 
entering the regi.on. In an actual blade channe.l, particles which en- 
ter the near blade region normally collide with, the adjacent blade 
surface and rebound back into the flow stream worth an increased vel- 
ocity lag and angular deviation relative to the gas. In practice, 

LDV electronics systems are generally designed Ito descriminate against 
such particles either because of their large siise or because their vel- 
ocities are significantly different from the me. ; an velocity of surrounding 
particles. 
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GAS FLOW RESULTS 

The results of numerically integrating eqn. (1) through the so- 
lution region of figure 2 are shown in figures 3 and 4. Figure 3 
presents the tangential variation of the gas velocity magnitude at 
several meridional locations. Tangential profiles were determined 
along every vertical grid line of figure 2; however, in order to min- 
imize confusion, only profiles at \ inch increments have been plotted. 
Since the tangential limits of each velocity profile on figure 3 corres- 
ponds to the limits of its grid reference line of figure 2, the two 
figures must be used in conjunction with each other to determine the 
meridional position of each profile. In general, the profiles progress 
from the upper left to the lower right as the meridional stream length in- 
creases. Each intrablade profile includes the blade surface velocity as 
its limit points. The results indicate that the velocity on the suction 
side is generally greater than the pressure side, at the same meridional 

position, and that the flow is generally well behaved except for local 
flow reversals near the trailing edge of the upper blade. The intrablade 
Mach numbers range from 0.125 - 0.891. 

Tangential profiles of the angle of the gas velocity vector are 
presented in figure 4 at the same meridional positions as figure 3. 

Flow angles are measured from the meridional direction. Flow approaching 
the blade channel is represented by the profiles on the upper right 
portion of figure 4, and the profiles progress downward and to the 
left as the meridional streamlength increases. It's seen that the gas 
leaves the flow region at an angle of approximately -62.5 degrees, and 
that the flow angle is relatively insensitive to the tangential position 
in the channel except near the leading and trailing edges of the blade 


channel. 
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PARTICLE FLOW RESULTS AND DISCUSSION 

The results of numerically integrating eqhs. (3) and (4) through 
the solution region are summarized in figures 5-13. Figure 5 pre- 
sents several trajectories of 1 micron diameter particles with a mass 
density of 1 gm/cc. The trajectories are typical of those obtained 
throughout the study and represent the approximate range of trajectories 
examined which did not intersect either the upper or lower near-blade 
region. The partial trajectories shown downstream from the blades 
are those of particles which pass through the blade channel directly 
above the one investigated. 

The influence of particle radius and mass density is shown in 
figures 6 and 7, respectively. An increase in either parameter causes 

an extension, or lagging, of the trajectory and this results In a 

corresponding increase in the velocity lag and angular deviation be- 
tween the particles and the gas. While the results shown are for a 
particle mass density and diameter of 1 gm/cc and 0.5 micron, re- 
spectively, the results are typical of all other cases examined. 

It f s seen that the particle path is more sensitive to an increase in 
radius than to a corresponding increase in mass density. This is 
justified by noting from eqns . (3) and (4) that if Stokes flow pre- 
vailed (C D = 24/Re), then p p would appear to the first power in the 
denominator of both right hand terms and the particle radius would 
appear to the second power. Under these conditions a doubling of the 

particle size would have the same effect on the trajectory as a 4 to 1 

increase in mass density. While Stokes flow does not always occur in 
this study, the above analogy remains essentially valid since the drag 
expression given by eqn. (5) is based upon the Stokes drag law. 
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The tangential variation of particle velocity ratio and particle 
angular deviation at several meridional positions is illustrated in 
figures 8-11 as a function of particle radius and mass density. The 
abscissa is the normalized tangential position and is 0 at the suction 
surface and 1 at the pressure surface. Particle velocity ratio is de- 
fined as the magnitude ratio of the particle velocity to gas velocity 
V/U, and particle angular deviation $ is defined as (j) = T^ ese 

parameters characterize the dynamic behavior of a particle at a par- 
ticular point in a fluid flow. Particle velocity lag is defined as 
ju-V |/U and is a related alternate parameter. Such a parameter does 
not, however, indicate whether the particles are moving faster or slower 
than the surrounding gas. LDV applications typically require that 
.99 £ V/U £ 1.01 and 4> <.5 degree. The tangential profiles of V/U and $ do 
not extend the full tangential height of the stream channel because of 
intersection of the outer trajectories with the near-blade region and the 
channeling of trajectories with increasing meridional streamlength and 
particle radius. The channeling effect is particularly evident in fig- 
ures 8 and 10. 

The results indicate that an increase in either particle radius 
or mass density causes an increase in velocity lag and angular ' 
deviation. As figures 6 and 7 suggested, particle tracking is more 
sensitive to par-tide radius than to mass density. Figures 8-11 in- 
dicate that doubling the particle radius causes approximately the same 
change in V/U and throughout the region as a 4 to 1 increase in the 
mass density. Therefore, the particle trajectory would remain essen- 
tially unchanged if the particle diameter was halved while at the same 
time increasing its mass density by a factor of 4. The figures 
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also indicate that tracking is better in the region downstream from 
the blade channel than in the channel itself. This trend results from 
the deceleration of the gas in the latter portion of the stream channel 
(fig. 3) and the gradual catching-up of the particles relative to the 
gas. When X=2.26 in., figs. 8 and '9 indicate that in the region adjacent 
to the suction surface the particles are actually moving faster than the 
gas flow (V/U>1). 

Figures 12 and 13 illustrate the tracking capability of three 
different size particles through the solution region. Tracking is 
graphically shown by lines of constant velocity ratio and angular 
deviation. The results indicate that the steepest gradients in V/U 
occur near the leading and trailing edges of the blade channel. The 
high leading edge gradient is caused by the sudden contraction of 
the blade channel in that region, and the high gradient just downstream 
from the blade channel is due to the sudden enlargement of flow area 
at the trailing edge. 

The regions of greatest velocity lag occur generally in several 
well defined regions, or pockets.’ Two pockets are adjacent to the 
leading edges of the blades in the region of greatest curvature. 

Another pocket lies approximately two thirds of the way through the 
blade channel and adjacent to the pressure surface. Particles in these 
regions are moving slower than the gas and are attempting to catch-up. 

The remaining regions of high velocity lag are near the trailing edge 
of the upper blade surface and in the narrow, high V/U gradient region 
immediately downstream from the blade channel. Because of the gas de- 
celeration in these regions, the particle velocity actually exceeds 
the gas velocity. 



Minimum velocity lag generally occurs near the entrance to the 
flow region where the gas-particle mixture was introduced, and in a 
region just prior to the trailing edge of the channel where the 
particles have caught-up to the decelerating gas. Another favorable 
region is seen to lie in the narrow high gradient region downstream 
from the channel. 

CONCLUSIONS 

The results presented above have provided information regarding 
the dynamic behavior of micron size particles in the two -dimensional 
flow field between blades of a stationary turbine cascade. The re- 
sults indicated that particle velocity lag and angular deviation is 
increased whenever particle size or mass density is increased, and 
that these tracking parameters are more sensitive to a change in 
particle radius than to a corresponding change in mass density. It 
was found that a 2 to 1 increase in particle diameter decreased 
particle trackability by approximately the same amount as a 4 to 1 
increase in particle mass density. Therefore, a 0.5 micron diameter 
particle with a mass density of 4 gm/cc behaved essentially the 
same as. a 1 micron particle with a mass density of 1 gm/cc. The 
greatest velocity lag occurs in several regions near the entrance 
and exit of the blade channel and the minimum lag occurs in the 
region of decelerating gas flow prior to the blade channel exit. 
Results indicate that LDV applications employing 1 gm/cc tracer 
particles with diameters greater than approximately 1 micron, or 
0.5 micron diameter particles with mass densities greater than 4 
gm/cc would experience velocity and angular deviations generally 
greater than 2 percent and 1 degree, respectively. 
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It should be emphasized that this study did not investigate par- 
ticle flow in the boundary layers of the blade channel nor did it 
examine particle collision phenomena with the blade surfaces. There- 
fore, it is conceivable that there may be particles within the blade 
channel which have previously collided with a blade surface and possess 
a velocity and flow angle . significantly different from the surrounding 
undeflected particles. Several cases were numerically investigated in 
which particles were introduced into the blade channel with velocities 
and flow angles significantly different from those of the surrounding 
gas. The results generally showed that particles below 1 micron in 
diameter and less than 4 gm/cc adjusted very rapidly to the gas flow 
field around them. The conclusion is that unless the particles are 
very large and dense, the dynamic behavior of blade impacted particles 
is essentially indistinguishable from the behavior of the surrounding 
undeviated particles. 
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FIGURE 1 Blade-To-Blade Solution Region 
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FIGURE 2 Blade-To-Blade Solution Region 
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FIGURE 3 Gas Flow Velocity Profiles Across Stream Channel, 

m = 6.68 x 10~ 3 lb/sec, p = 7.403 x 10 lb/ft 3 , 
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T = 523.7 R, y - 1-4 
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FIGURE 4 Gas Flow angle Profiles Across Stream Channel, 

m = 6.68 x 10" 3 lb/sec, p = 7.403 x 10 -2 lb/ft 3 , 
T = 523.7 R, y = 1-4 
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FIGURE 5 Particle Trajectories Through the Solution Region 
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FIGURE 6 Particle Trajectories as a Function of Particle Radius 
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FIGURE 7 Particle Trajectories as a Function of Particle Mass Density 
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FIGURE 12a Contours of Constant Parti cle-To-Gas Velocity Ratio 


p =1 gm/cc 
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FIGURE 12b Contours of Constant Particle-To-Gas Velocity Ratio 
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FIGURE 13a Contours of Constant Angular Deviation, 

P = 1 gm/cc 
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Contours of Constant Angular Deviation. 
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FIGURE 13c Contours of Constant Angular Deviation, 


p =1 gin/cc 
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FIGURE 2 Blade-To-Blade Solution Region 
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FIGURE 6 Particle; Trajectories as a Function of Particle Radius 
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FIGURE 7 Particle Trajectories as a Function of Particle Mass Density 
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FIGURE 12a Contours of Constant Parti cle-To-Gas Velocity Ratio, 


p =1 gm/cc 
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FIGURE 12b Contours of Constant Particle-To-Gas Velocity Ratio, 

P p = I gm/cc 
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FIGURE 12c Contours of Constant Particle-To-Gas Velocity Ratio, 


P_ = 1 gm/cc 
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FIGURE 13a Contours of Constant Angular Deviation, 

P p = 1 gm/cc 
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FIGURE 13c Contours of Constant Angular Deviation, 


p =1 gm/cc 
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